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We examine random matrix systems driven by an external field in view of optimal control theory 
(OCT). By numerically solving OCT equations, we can show that there exists a smooth transition 
between two states called "moving bases" which are dynamically related to initial and final states. 
In our previous work [J. Phys. Soc. Jpn. 73 (2004) 3215-3216; Adv. Chem. Phys. 130A (2005) 
435-458], they were assumed to be orthogonal, but in this paper, we introduce orthogonal moving 
bases. We can construct a Rabi-oscillation like representation of a wavpacket using such moving 
bases, and derive an analytic optimal field as a solution of the OCT equations. We also numerically 
show that the newly obtained optimal field outperforms the previous one. 



PACS numbers: 05.45.Mt, 02.30.Yy, 03.65.Sq 



I. INTRODUCTION 

Controlling atomic and molecular processes by laser 
fields is one of current topics in physics and chemistry 
There are various control schemes applied to such 
processes: 7r-pulses Q .nonadiabatic transitions [H, adi- 
abatic rapid passage |4j .Stimulated Raman Adiabatic 
Passage (STIRAPJjl, H, 0] , pulse-timing control 0], and 
coherent control [3], etc. These strategies are known to 
work when the system to be controlled is rather simple or 
small. However, the system can be complex (Io| when we 
deal with highly excited states in large molecules or meso- 
scopic devices driven by electro-magnetic fields. Such a 
"complex" system in the limit of strong chaos is modelled 
by a random matrix Hamiltonian with a time-dependent 
external field [ll|, and the dynamics is well represented 
by multi-level-multi-level transitions with random inter- 
actions among energy levels. Although there are many 
works on the statistical properties fl2l Il3l. [Til Il5l . [lt| as 
well as the semiclassical properties [13, Qjl of eigen- 
values under the variation of an external parameter, few 
works have been published on the dynamical properties 
of such systems except several studies on nonadiabatic 
processes [13, [H|. 

Even for such complex systems, there exist mathemat- 
ical results showing complete controllability [2^, HH of 
general quantum systems with discrete spectrum under 
certain conditions. The existence of an optimal field 
is proved by optimal control theory (OCT) [24 ]. which 
is a powerful tool to obtain an optimal field and has 
been studied for various dynamical systems [25L l26| . For 
the purpose of steering quantum states, many numer- 
ical schemes with monotonically convergent algorithms 
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0, H3, HI have been developed based on OCT. In gen- 
eral, OCT for quantum states provides sets of nonlinear 
differential equations (OCT equations) which are solved 
by iterative procedures. For complex systems with many 
degrees of freedom, however, the optimal field often be- 
comes too complicated to analyze the dynamical pro- 
cesses involved. In addition, the computational cost be- 
comes significantly heavy when we apply OCT to real- 
istic problems with many degrees of freedom. Analytic 
approaches can be a good strategy to complement this 
annoying situation. 

One such analytic method for multi-level control prob- 
lems is STIRAP @, H, 13- 

Though it can accomplish 
perfect control, it assumes an intermediate state coupled 
to initial and target states, and uses a pair of external 
fields with slowly varying amplitudes. Rece ntly , we have 
proposed another analytic optimal field 29, 30] which in- 
duces a "direct" transition between random vectors in a 
random matrix system. The key idea of this approach is 
to describe the optimally controlled dynamics as a Rabi- 
like oscillation [2j], and our optimal field can be inter- 
preted as a generalized 7r-pulse |29j |. Though the deriva- 
tion and applicability of our analytic optimal field have 
been detailed in |3Q ] , there exists deficiency in our previ- 
ous formula because of several (unnecessary) assumptions 
for simplification. In this paper, we rederive an analytic 
optimal field with less numbers of assumptions and reex- 
amine its applicability to random matrix systems. 

This paper is organized as follows. In Sec. [Til we nu- 
merically investigate the multi-state control problem by 
OCT to show that, in some cases, the optimal field in- 
duces a smooth transition. According to this observation, 
in section IIII1 we introduce a Rabi-like representation of 
the controlled state with some modifications compared to 
our previous result [29L l30l ] . Employing this representa- 
tion, wc obtain a new analytic expression of the optimal 
field. In Sec. IIII El we confirm the applicability of the 
analytic field through the numerical integration of the 
Schrodinger's equation for random matrix systems. Fi- 
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nally, in Sec. lIVi we summarize this paper and give some 
discussions on the control problem of quantum chaos sys- 
tems. We mention some technical details in Appendix. 



II. OPTIMAL CONTROL IN RANDOM 
MATRIX SYSTEMS 

We present numerical results of controlled dynamics 
driven by an optimal field to see what kinds of dynamics 
are involved in random matrix systems. The random 
matrix Hamiltonian driven by a time-dependent external 
field e(t) is written as 



H[e(t)] =H + e(t)V, 



(1) 



where Hq and V are random matrices subject to a certain 
universality class [ll[, i.e. Gaussian Orthogonal Ensem- 
ble (GOE), Gaussian Unitary Ensemble (GUE), etc. 

It is well known that a strongly chaotic system does 
not have any constant of motion except the total en- 
ergy flfj| . where the typical quantum states are ran- 
dom vectors. Thus, it is appropriate to choose initial 
and target states as random vectors. If we choose a 
certain ortho-normalized basis, a random vector in TV- 
dimensional Hilbert space is represented by a set of ran- 
dom complex numbers {cj}. If such a vector has neither 
special symmetry nor correlation, only the constraint im- 
posed is the normalization condition, 



JV 

Ei 



(2) 



Then, the normalized probability density for a variable 
y = | cj | 2 is given by 



Pn(v) dy = Nexp(-Ny) dy, 



(3) 



when N is sufficiently large 

The actual procedure to numerically obtain an optimal 
field is as follows: The Hamiltonian {T]) is constructed by 
generating two random matrices, Hq and V, with N x 
N elements, where the scales of them are determined 
so that the averaged eigenvalue-spacing AE of Hq and 
the variance AV of off-diagonal elements of V are both 
unity. Next, we define an initial state |$o) aid a target 
state |$t) as random vectors satisfying the distribution 
([3|). Then, for a fixed target time T, the optimal field 
e(t) is obtained by solving the OCT equations which are 
detailed in Sec. Ill Al 



A. Zhu-Botina-Rabits Scheme of OCT 

There are many effective methods to solve OCT equa- 
tions for quantum systems [lj. In this section, we use 
a method introduced by Zhu, Botina, and Rabitz [27| 
(ZBR-OCT). Our goal is to determine the optimal exter- 
nal field e(t) by which a given initial state |$o) is steered 
to a given target state |3>t) at a target time T. According 
to ZBR-OCT, we introduce a functional J(e(t), \<fr(t))) 



J(e(t), |0(t)» = Jo - a [ [e(t)} 2 dt - 2Re (</>(T)|<f T ) f 
Jo Jo 



("<T>|<i>,> l' T (x(t)\§- t -^§^Mt))dt 



(4) 



r 



where T and a are given parameters representing the tar- 
get time and the penalty factor, respectively. The quan- 
tum state \4>(t)} satisfies the initial condition, |</>(0)) = 
|$o)- The first term in the right-hand side is the final 
overlap, 



Jo = |^(T)I*t)| S 



(•5) 



The second term is the penalty term which minimizes the 
amplitude of the optimal field. In the third term, a La- 
grange multiplier |x(i)) is introduced to give a constraint 
that \4>(t)) satisfies Schrodinger's equation, 



in-\<t>(t))=H[<t)]\m- 



(6) 



On the other hand, Schrodinger's equation for |x(*)) is 



ih-\ X (t))=H[e(t)]\x(t)), 



(7) 



and the boundary condition |x(T)) = are obtained 
by "differentiating" the functional with respect to | </>(£)) 
and \(f>(T)). For the Hamiltonian ((U), the variation of J 
with respect to e(t) gives an expression for the optimal 
field, 



e(t) = —lm[m)\x(t))(x(t)\V\m] 
an 



(8) 



This is a self-consistent expression for the optimal field, 
and to obtain its actual value, we have to simultaneously 
solve the nonlinear coupled equations, ([6]), (0) and ([8]) 
(OCT equations). ZBR-OCT is one such method which 
numerically solves the OCT equations with iterative pro- 
cedures 12711. 
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FIG. 1: Numerical results of controlled dynamics between random vectors in a 128 x 128 GOE random matrix system. The 
optimal field is obtained through the iterative procedure given by Zhu, Botina, and Rabitz |27| ] for the target time T — 10 and 
the penalty factor a = 1. (a) Convergence property of functional values, Jo (solid curve) and J (dashed curve), as a function of 
the iteration step, (b) Optimal external field e(t). (c) Power spectrum of s(t). (d) Time evolution of the overlap \(&T\tp{t)}\ 2 
(the magnified curve near t = T is shown in the inset). 




FIG. 2: The same as Fig. [T] except that the target time is longer (T = 100) and the penalty factor is larger (a = 10). 
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FIG. 3: Time evolution of the overlaps {P x {t)) (solid curve), (Pj,(i)) (dashed curve), and (P^+ X (t)) (thick curve) defined by 
Eqs. (TlT}, (12]) and (14]). ( a ) and ( b ) correspond to Fig. [TJ (T = 10 and a = 1) and Fig.[2](T = 100 and a = 10), respectively. 



B. Numerical Results 



We show controlled dynamics driven by numerically 
obtained optimal fields for a 128 x 128 GOE random 
matrix Hamiltonian. The quantum state \ip(t)) with 
the initial condition \ip(0)) = |$q) evolves according to 
Schrodinger's equation ©. We have chosen Hq and V so 
that AE = AV = 1. In other words, the energy values 
are shown in unit of AE, and the unit of time is h/ AE. 
Then, the field strength e(t) is shown in unit of AE/ AV. 

In Fig. [TJ we show the result with the parameters 
T = 10 and a = 1. The target time T = 10 is com- 
parable to the minimum time To rs 2-7T (h/AE) which is 
necessary to resolve each energy level from its adjacent 
levels. In Fig.QJa), the functional values Jo © and J 
are shown as the solid and dashed curves, respectively. 
They appear to converge after several ten steps. The fi- 
nal overlap Jo is 0.89 after 100 iterations. The optimal 
field e(t) is shown in Fig.[TJb) as well as its Fourier spec- 
trum in Fig. [TJc) . Figure [TJd) shows the time evolution 
of the overlap | ($^1^(0) | 2 with its magnification near 
the target time in the inset. 

Figure [2] shows the result obtained for the parameters 
T = 100 and a = 10, which is the case of a relatively long 
target time compared to to. The values of Jo and J, the 
optimal field e(t), its Fourier spectrum, and the overlap 
|( ( I > t|V ; (0)| 2 are shown as in Fig. [TJ In this calculation, 
the final overlap Jo is 0.93 after 100 iterations. 

In both cases, the overlap |($t , |V'(*))| 2 as a function 
of time t remains small until t is close to the target time 
T . In multi-state quantum dynamics, even if no external 
field is applied, an auto-correlation function (ip(0)\ip(t)} 
can rapidly decay by dephasing among dynamical phases 
of Hq. This is the reason why the overlaps ( ( 1>t|V ; (0) 
in Fig. [TJd) and Fig. [^d) rapidly grow up to the final 
values near t = T. In other words, ($r|^(t)) decays 
quickly when t deviates from T. 



C. Observation of Smooth Transitions 

Since we want to concentrate on transitions induced 
by e(t) only, it is necessary to remove the contribution 
from dephasing by Hq. This is nothing but the proce- 
dure of the interaction picture in quantum mechanics 
[3lj |. We define time-dependent quantum states related 
to |$o) and |$ T ) by 

\M*)) = #b(t, 0)|*o>, Ixo W) = Ua(t,T)\*T), (9) 

where Uo(t2,ti) represents a propagator from t = t% to 
t = ti with respect to the unperturbed Hamiltonian Hq. 
We call these states "moving bases." In the following, we 
analyze the optimally controlled dynamics through these 
time-dependent states. 

If we introduce projection operators associated with 
these states by 

P^t) = \Mt))(Mt)l P x (t) = \ X o(t))(xo(t)\, (io) 

the probabilities such that \ip(t)) is found in these states 
are written as 

(P^t)) EE Mt)|J^(t)>, (11) 

(p x (t)) = m)\Px\m)' (12) 

These values are more appropriate quantities to observe 
the multi-level-multi-level transition dynamics compared 
to the bare overlap | (*&t" | "0 (*) ) | 2 as shown below. In ad- 
dition, we introduce another projection operator 

p (t) = hit) + p x (t) ~ p^{t)p x {t) - p x (i)p^t) 

l-tr[P x (t)P,(t)} 

(13) 

which represents projection onto a subspace defined by 
a linear superposition of \<fio(t)) and \xo{t))- We can 
prove that this is a projection operator by using 1 — 
tt[Pi(t)P x (t)] = 1 - \(Mxo)\ 2 and P+(t)P x (t)p4(t) = 
l(0o|xo)| 2 ^(O> etc Then, the quantity, 

(iW*)) = W*)I*W*)M*)>. ( 14 ) 
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represents the probability that the quantum state \ijj(t)) 
is found on the subspace. 

In Fig. [31(a) and (b), we show the overlaps (probabil- 
ities), (P^(t)), (P x (t)), and (P^ +X (t)), calculated from 
the results in Figs. Q] and All the curves in Fig. [3^b) 
are smoother than those in Fig. [3fa). It is also worth 
noting that (P^+ X (t)) stays close to unity for all the time 
in Fig. [lb). 

From other ZBR-OCT calculations for random matrix 
systems, we found that the ZBR optimal field induces 
a transition from |$ ) to |$t) nearly within a subspace 
spanned by \(j>o(t)) and \xo(t)) when the target time T is 
sufficiently large. Based on this finding, we will develop 
an analytic approach for the optimal field in the next 
section. 



III. ANALYTIC APPROACH FOR 
CONTROLLED DYNAMICS 

The Rabi oscillation in a two-level system has been 
studied in detail [2j. According to such previous works, 
we can represent a wavefunction as a linear combination 
of two eigenstates \tp\) and \tp2) with eigen-energies E\ 
and E2, 

m)) = A{t)\<p x )e E ^ in + B(t)\^ 2 )e E ^ lh . (15) 

Here the coefficients A(t) and B(t) are slowly oscillating 
functions with a Rabi frequency under the rotating-wave 
approximation (RWA) Q. 

In this section, we show, under certain conditions, that 
a Rabi-like description becomes valid even for multi-level 
quantum systems where the wavefunction is described by 
the time-dependent states in Eq. (jjJJ). We call these states 
"moving bases" instead of eigenstates. This is equivalent 
to considering the case where the controlled state remains 
in the subspace spaned by these moving bases over a 
whole period of the dynamics. With the help of OCT, we 
conversely obtain an analytical expression for the optimal 
field to induce the smooth transition we found in Fig. [3] 
for the multi- level dynamics. 



A. Rabi-like Representation 

In the previous section, we have defined the moving 
bases |</>o(i)) an d lxo(*)) © i n order to observe smooth 
transitions in OCT calculations. Unlike the two-level 
case, however, these states are not always orthogonal to 
each other. The inner-product between them is written 
in the form, 

(Mt)\Xo(t)} = ($ \U(0,T)\<S> T ) = ze^sinO, (16) 

with < 9 < vr/2 and < 9 < 2n. Note that 
the inner-product does not depend on t, i.e. con- 
stant in time, despite that the moving bases (|0o(t)) and 
\Xo(t))) rapidly change their "directions" according to 
the Schrodinger's equations. 

In our previous works [29l . [30j , we have used an assump- 
tion that O = 0. Actually we can remove this assumption 
by introducing an orthogonal pair of the moving bases as 
(Schmidt decomposition) 

\Mt)) = |0o(*)>, (17) 
I - , A . |xo(t))-»e i9 sin9|fo(*)) (m , 

IXOW) - K ■ (J-o) 

cos O 

These are our new moving bases which are orthogonal to 
each other and will be used below. 

We introduce a Rabi-like representation of the quan- 
tum state \4>{t)} with an initial condition |^>(0)) = |$o) 
by a linear combination of the new moving bases \(f>o(t)} 
and \xo(t)), 



^(t)) = A(t)\Mt)) + B(t)\xo(t))- 



(19) 



The coefficients, A(t) and B(t), must satisfy a normal- 
ization condition, 



\A(t)\ 2 + \B(t)\ 2 = l. 



(20) 



If this representation is valid, A(t) and B(t) satisfy the 
following differential equations, 



> - 


ih 


(Mt)\v\Mt))A(t) - 


-(Mt)\V\xo(t))B(t) 


> - 


£(t) 

ih 


'(Ut)\v\Mt))A(t)- 


H (xo(t)\V\xo(t))B(t) 



(21) 
(22) 



r 



B. Rotating Wave Approximation 



The rotating-wave approximation (RWA) means drop- 
ping rapidly oscillating terms in differential equations [2| . 
Since RWA is applicable to non-degenerate multi-level 



systems 0, we will apply RWA to Eqs. (JUD and ^ to 
solve them approximately. 

To justify the use of RWA, we introduce three integrals, 

In(t) = ~ / e(t)(Mf)\V\xo(f))df, (23) 
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FIG. 4: Time evolution of \In(t)\ (thick line), \I<f,(t)\ (solid line), and |/ x (t)| (dashed line), defined by Eqs. lf23]l, ((24)l , and (|25]| , 
respectively, (a) and (b) correspond to Fig. [T](T = 10 and a = 1) and Fig. [2](T = 100 and a = 10), respectively. 



= ^ ^(^(Ol^oC*')}*', (24) 
= T A(*)<Xo(OMxo(*'))d*'- (25) 



Note that the generalized pulse area [29| is represented by 
2|Jjn(T')|. The numerical results in Figs. [Hand [2] are used 
to calculate these integrals explicitly, and their absolute 
values are plotted in Fig. 2J Since /</>(£) and I x {t) are 
almost zero, it is appropriate to assume that 



I (j> (t)=I x (t) = 0, 



(26) 



which corresponds to RWA. Furthermore, for a long tar- 
get time, we can employ the following form 



ih(t) = nt, 



(27) 



which means that, with use of RWA, the integrand in 
Eq. (|23|) is nearly constant in time. 



Under these approximations (f^Bj) and l|27p. the differ- 
ential equations, (j2l"j) and (f22)) . are simplified to 



—A(T) = -iVLBti), —BIT) = -in* Ait). (28) 
at at 



For the initial value problem with the conditions, A(0) 
1 and B(Q) — 0, we obtain a solution, 



A(t) = cos[|n|t], £(t) = -^sin[|0|t]. (29) 



Defining the phase of 17 as 

W\" 

we finally obtain 



(30) 



\m = \Mt)) cos[\m - «^ixo(i)> M\m = coa[l ^ e] \Mt)) - te ~^T |t] Mt))- (31) 

This state oscillates with the frequency |f2| between \<fio(t)) and |xo(^))j as well as between \(j>o(t)) and |xo(i))- 



C. Analytic Field 

According to the ZBR-OCT scheme in Sec. Ill Al the optimal external field is represented by Eq. ([8]) given the 
forward evolving state \4>(t)) and the backward evolving state |%(t)) are prepared with the boundary conditions, 

10(0)) = |$ ), IxCO) - |$t>. (32) 

Under the approximations (|26|) and (|27|) . \<fi(t)) has been already given in Eq. (|3Tj) . and |x(*)) is written as 

| x (i» = -ie"|A,(«))sta[|!l|(«-T)-e] + |xo(t))cos[|!!|(t-T)-e] (33) 

cos B cos fc) 
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FIG. 5: Numerical results (squares and circles) by OCT with the penalty factor a = 1 are compared to our theoretical values: 
(a) Final overlap Jo. (b) Averaged amplitude s of the external field. The solid (dashed) curve corresponds to the case of 
sine = t/tt/4/V (8 = 0). 
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FIG. 6: The same as Fig. [5] except that the penalty factor is larger (a = 10). 



The overlap between these states is 

(ct>(t)\x(t))=ie ie MMT + e}. 



(35) 



Substituting Eqs. $tt$ and |(33J) into \<j>(t)) and |x(£)) in 
Eq. (T5J) , we obtain the external field, 



2afi 



e i9 (Ut)\V\Mt)) • (36) 



By calculating Iq (T) for the external field ([36]) and using 
Eq. (|27|) at t = T, we obtain an equation for O, 

n= Bin[2(|n|r + e)] [e ^ 2 + e -^ 2] (37) 

where V 2 and W 2 are defined by averages of transition 
elements, 



V 2 



w z 



1 

f 



T 



(Mt)\v\m) 



(Mtwim) 



dt, 



dt. 



(38) 
(39) 



Note that \W 2 \ becomes small compared to V 2 when the 
system is sufficiently large without special symmetry (see 



Eqs. (jA20l) and ([All]) ). Then, we obtain an equation for 
|f2| with use of Eq. ([50]). 



ini = 



v 2 sin[2(|n|r + e)] 

4ah 2 



(40) 



The solutions of this equation are obtained from the 
crossing points between y = x/T and y — (K/2) sin[2(a;+ 
8)] where x = \fl\T and K = V 2 /(2ah 2 ). 

To illustrate the effectiveness of our result, we calculate 
the final overlap Jq from Eq. (|35| 



Jo = Mt)\x(t))\ 2 =vAMT + G] (41) 

and the average amplitude of the optimal field from 
Eq. 



(42) 



These estimates match well with the numerical results 
(squares and circles) as shown in Figs. [5] and [5J Since 
the inner product between N dimensional random com- 
plex vectors is ~ yJir/AN as an average (see Eq. (|A19[) ) . 
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we have used \(4>o(t)\xo(t))\ = sin 6 ~ yJir/AN (solid 
lines) . For comparison, we also show the results for = 
(dashed lines), which correspond to our previous results 
|29l l30j ] . It is obvious that the new analytic field outper- 
forms the previous one. 



D. Perfect Control 



In the ZBR-OCT scheme (Sec. Ill A|) . we need a finite 
value of the penalty factor a in order to avoid numerical 
unstability. In the analytic approach, however, the limit 
a — > can be taken safely in Eq. (|40|) . In this case, the 
solution of Eq. (|4"0")) becomes 



101 -» fL 



1 /7T 

T \2 
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(43) 



where m is an integer. It is easily shown that this gives 
J = 1 horn Eq. (|H|). From Eqs. J3B]| and (gOJ), the 
external field is obtained as 



E (t) 



-Re 



(Xo(*)|^o(*)> 



(44) 



Since this equation does not contain a, it is different from 
other non-iterative optimal fields (32j. It can be shown 
that this field actually achieves perfect control in the 
limit of T, N — > 00. We give the proof in Appendix I A 41 



E. Application of the Analytic Field 

Although it is theoretically exact for T, N — > 00, our 
approach is applicable to the cases with finite T and N 
as already given in Figs. and [B] 

In Fig. [71 we show time evolutions of the overlaps with 
T = 10 (a) and T = 100 (b). The initial and target 
states are the same as in Fig. [3] Unlike the results by 
OCT (Fig. [3]), the probability on the subspace spanned 
by |0o (*)) and |xo(*)) decreases monotonically. This is 
because the analytic field ([44|) was obtained under the 
conditions T, N — > 00 while the numerical calculations 
were performed for finite T and N . 

The performance of the analytic optimal field can be 
easily seen by plotting the final overlap Jo for various 
values of N and T (Fig. [8]). The errorbars in this fig- 
ure represent the normal deviation of Jo obtained from 
calculations for 100 different samples of the Hamiltonian 
and state vector. By comparing with our previous re- 
sult, Fig.l in [2!^ ]. our new analytic field outperforms the 
previous one for the intermediate T and N as expected. 
This is a nice feature when we consider the application 
of this method to non-limiting cases. 



From the analysis of OCT calculations, we showed that 
optimally controlled states remain on a subspace spanned 
by two "moving bases" when the target time and sys- 
tem size are both sufficiently large. According to this 
observation, we developed a new method to solve OCT 
equations and to obtain an analytic expression for the 
optimal field. Finally, it was numerically shown that the 
analytic field actually steers the quantum states in ran- 
dom matrix systems. The difference from our previous 
result is that we have taken new moving bases which are 
exactly orthogonal, and the newly obtained analytic field 
outperforms the previous one for the intermediate target 
time and system size. 

Our analytic field (|44|) is a generalized 7r-pulse [2^] in 
multi-level systems, which is realized because of certain 
randomness in the elements of the Hamiltonian and state 
vectors. The amplitude of the pulse becomes smaller 
when the target time T is larger since an effective pulse 
area should be a constant tt [29| . Although our controlled 
dynamics seems to be antithetical to the molecular pro- 
cesses induced by intense laser fields [H, |34|, this does 
not necessarily mean that our approach is not applicable 
to those systems. This is because such dynamics driven 
by the intense laser field can be included in the "unper- 
turbed" Hamiltonian Hq . If the system becomes strongly 
chaotic by the laser field, such a situation is even prefer- 
able for the prerequisite of our approach using random 
matrix Hamiltonians. 

The quantum targeting problem studied in this paper 
was solved analytically for random matrix systems, while 
it is known that classical targeting problems [HI, [3(| are 
difficult to be solved for strongly chaotic cases. This is 
because there is sensitivity of trajectories with respect 
to initial values. Our result for quantum systems thus 
seems to break the naive quantum-classical correspon- 
dence. It is important to clarify how the correspondence 
is recovered in the semiclassical limit [13, HH . 

Recently Gong and Brumer showed that coherent con- 
trol works for a quantized kicked rotor [39l . |40J , a typical 
"quantum chaos" system, whereas our concern was op- 
timal control of quantum states in random matrix sys- 
tems. Optimal control for quantum chaos systems, es- 
pecially weakly chaotic systems, is another interesting 
subject which should be pursued [30ll4l|. 
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FIG. 7: The same as Fig. [3] except that we have used the analytic optimal field Eq. (144 \ . (a) and (b) correspond to T = 10 
and T = 100, respectively. 
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FIG. 8: The final overlap Jo obtained by the analytic optimal 
field. We show the results for various system sizes according 
to the following procedure: At first, we generate two random 
matrices (Ho and V) and two random vectors (|3>o) and |$t)) 
for a system size N . If we choose a target time T, the analytic 
field is given by Eq. (|44]l . The quantum state \ip(t)) at t = T 
is obtained by numerical integration of Schrodinger's equation 
under the field with the initial condition ^(0)) = |$o). Fi- 
nally, the final overlap Jo is calculated by \(&t\tP(T)}\ 2 . Each 
curve in the figure is obtained as an ensemble average over 100 
different realizations of random numbers. Compare this with 
Fig. 1 in the previous study [29(] (only the data of N = 32 is 
shown as a dotted curve in this figure for comparison). 



APPENDIX A: EIGENSTATE 
REPRESENTATION 



1. Preliminary 

Though our results in the main text do not depend on a 
particular representation, in this appendix, we describe 
the controlled dynamics and the analytic optimal field 
Eq. (|4"4"|) by using the eigenstate representation of Ho, 
and prove that perfect control is achieved by Eq. (|4"4")) . 

We introduce eigenstates \<fj) of H corresponding to 
eigenvalues Ej where {|<^j)} constitute an orthonormal 



JV 



N 



l $ o) = ^Cj\<p S ) 
3 = 1 



|*r) = X; d il^) ( A1 ) 

3=1 



are represented by random complex numbers {cj} and 
{dj} satisfying normalization conditions, 



N N 
3=1 3=1 



(A2) 



The matrix elements of V are defined by using { | tfj)} 

V jk = (<Pj\V\<p k ). (A3) 

These quantities, {cj}, {dj} and {Vjk}, are assumed to 
be uncorrelated among them. 

2. The Analytic Optimal Field 

The moving bases, Eq. ([5]), satisfying |^(0)) = |$o) 
and \x(T)) = |$r) can be written as 



JV 



3 = 1 
JV 

3=1 



(A4) 
(A5) 



where T represents the target time. In general, these 
states are not orthogonal to each other as shown in 
Eq. (fl6|) . and the orthogonal (new) moving bases, 
Eqs. (fTT)) and ([T5jl . are constructed as 



\Mt)) = Y. c ^ Eit/in 

3 = 1 
N 



(A6) 



dje -E s T/in_ ie ie CjShie ^^ E _ t/in 



\ipj)e 



3=1 
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JV 



3=1 



(A7) 



Substituting these states into Eq. ([44]) . we obtain the 
eigenstate representation of the analytic optimal field, 



2hQ m 



N N 



^ ^ ^ ^ feCfc exp \ 

7=1 fe=l ' 



(A8) 



3. Sum of Random Variables 



For large N, we obtain 



JV 



J2\V kj \ 2 \ Cj \ 2 ~ JVA4(|V fci | a | Cj | 2 ) 

3=1 



E^l 2 !^ 



= X(|^-| 2 ) + 0(l/ViV), (A14) 



iVX(|4| 2 |^| 2 ) 



fc=l 



- M(\V kj \ 2 ) + 0(l/VN), (A15) 
where we have used a basic relation 



Suppose a probability variable F is defined by a sum 
of independent probability variables Xj (j — 1, . . . , n), 



y = Xi + x 2 



A- 



(A9) 



with an expectation A4(Xj) = A4x (> 0) and a vari- 
ance a 2 (Xj) = ax 1 - When n is sufficiently large, it is 
known from the central limit theorem that Y is normally 
distributed. The expectation and variance are 



M(Y) = nMx, <J 2 {Y) = na x ' 



(A10) 



Then, the expectation nMx can be used as an approxi- 
mate value of Y since the relative standard deviation 



Orel 



M{Y) MxV^ 



0(1/VH) (All) 



vanishes for n — * oo. 

Using the above basic knowledge, we estimate approx- 
imated values for sums of V k j , Cj , and d k in the following. 
The coefficients Cj and dj are independent random num- 
bers subject to the distribution function ([3]), i.e., 



M(\ Cj \ 2 ) = M(\dj\ 2 ) = -^r, 
a 2 (| C ,| 2 ) = a 2 (\d j \ 2 ) = ± 



(A12) 
(A13) 



^(Xi^a-.-X,,) =M(Xi)M(X 2 )---M(X n ) (A16) 

for independent probability variables Xj. 

Applying the central limit theorem to a sum of complex 
variables Z n = z\ + Z2 + - ■ - + z n = X n +iY n with M(zj) = 
and M(\zj\ 2 ) = a 2 , we have 



P(X n )P(Y n ) dX n dY r 

1 



■ exp 



x n 2 + Y n 2 



dX n dY n . (A17) 



Then, the average magnitude of \Z n \ can be calculated 

as 



M(\Z n \) = J \Z n \P(X n )P(Y n ) dX n dY n = 

(AIM 

From this relation, the inner product, Eq. (|16p . is esti- 
mated as 




(A19) 



In the same manner, the average transition elements 
and ([55)1 are estimated as 



E E k V ^k * M{\V kj \ 2 ) + 0(1/N), 



3 kjtj 



W 2 = ( E c i^i<*> I + 2 E E C K- l^l 2 did* - 0(l/iV), 

3 k<j 



(A20) 
(A21) 



in the limit T —> oo. Thus, we can ignore W 2 for T, N 
oo. 



4. Controlled State 

We give a proof that a quantum state driven by the 
analytic field, Eq. (|44f or (|A8|) . actually shows a smooth 
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transition between the initial and target states. We as- 
sume that the size N of the random matrix Hamiltonian 
H and the target time T are both large enough. 

To see the dynamics induced by the field, we represent 
a quantum state in the eigenstate representation, 

ltKt)>=E *(*)lw> eB * t/W ' (A22) 

3 

which satisfies the initial condition |^(0)) = |$o)- From 
the Schrodinger's equation driven by the optimal field 
(|A8|) . we obtain the following differential equations for 
dj(t), 

ak(t) = -iJ2 n kjaj{t), (A23) 

where VL k j is defined by 

n kj = ^^(c*d k +d*c k ), (A24) 



and we have used the rotating-wave approximation. We 
write a k (t) in the following form, 



a k (t) = A k cos fl m t + B k sin H m t, (A25) 



and A k and B k are determined as 



A k = c k , B k = -iJ2 77^" c i' (A26) 



from the initial conditions a k (0) = c k and (|A23() . 
Using the relations (|A14|) and (|A15|1 . we obtain 



j=jtk j^k j^k 



v-flkj,- j \- \Vkj\ 2 c*dj ^ \Vkj\ 2 \dk\ 2 ,nrt/An 

l-* j = k ^ — Cfc ^ ^ — ~ Cfe 0(1/iV) - (A28) 



Substituting (|A27jl into (|A25|1 . we obtain 

a k (t) = Cfc cos tt m t — idk sin tt m t, ( A29) 

and the right-hand side of the differential equation (|A23|) 
becomes 

— i £l k jO,j (t) = — O m ^Cfc sin f2 m £ + id k cos f2 TO i^ 

(A30) 

with use of (|A27|) and (|A28j) . Since Eq. (|A30|) is exactly 
the same as d k (t), we have confirmed that (|A29[) is the 



solution for the Schrodinger's equation driven by the an- 
alytical optimal field. 

The final expression (|A29[) shows that each a k (t) 
smoothly changes its value from c k at t — to d k at 
t = T as expected, and the overlap between \if>(t)) and 
|Xo(*)) is easily calculated as 

(Xo(t)\ip(t)) =-ismQ m t. (A31) 

This shows | (^tI^C^)) | — 1, i.e. perfect control is ac- 
complished at the target time T. 
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